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As a half-metal is spin-polarized at its Fermi level by definition, it was conventionally thought to 
have little proximity effect to an s-wave superconductor. Here we show that with interface spin-orbit 
coupling p x + ip y superconductivity without spin degeneracy is induced on the half-metal, and we 
give an estimate of its bulk energy gap. Therefore a single-band half-metal can give us a topological 
superconductor with a single chiral Majorana edge state. Our band calculation shows that two 
atomic layers of VTe or Cr02 is a single-band half-metal for a wide range (~0.1eV) of Fermi energy 
and thus is a suitable candidate material. 



Introduction: Possibility of Majorana fermions arising 
out of condensed matter system has aroused great inter- 
est in recent year*4 One class of systems where Majo- 
rana fermions can appear is the two-dimensional (2D) 
chiral superconductor which has a full pairing gap in the 
bulk, and Af gapless chiral one-dimensional (ID) states, 
which consists of Majorana fermions^l, at the edge. In a 
Af = 1 chiral topological superconductor (TSC), a single 
Majorana zero mode is bound to a vortex core^"^, giving 
rise to non-Abelian statistics which can be potentially 
useful for topological quantum computatiorP. The most 
straightforward way to realize such a chiral TSC is the in- 
trinsic p x + ip y superconductivity in spinless fermions^. 
The strongest candidate material for this superconduc- 
tivity, albeit a spinful version, is Sr 2 Ru04 7 , but the ex- 
perimental situation is not definitive^. Recently, there 
has been alternative proposals involving inducing s-wave 
superconductivity in material with strong spin-orbit cou- 
pling through proximit y e ffeclP^. It was pointed out 
in one of the proposals^ that spin polarized p x + ip y 
superconductor can be obtained in a ferromagnetic film 
through proximity to a superconductor. In this Letter, 
we further develop this line of approach and demonstrate 
through explicit calculations the feasibility of creating 
Majorana fermions in a half metal / conventional s-wave 
superconductor heterostructure. 

We consider the pair formation on a half-metal (HM) 
that is in proximity contact to an s-wave superconductor 
(SC). A HM, by definition, is spin-polarized at the Fermi 
surface 1 -, i.e. it is a metal for the majority-spin and an 
insulator for the minority-spin. Our proposal has two 
major advantages over other current proposals. Firstly, 
the Af = 1 TSC phase in our approach exists in a wide 
range of Fermi level. Secondly, we expect the proximity 
effect between SC and HM to be more robust than that 
between SC and semiconductor, due to better Fermi sur- 
face matching. It has been known that at the normal 
metal to s-wave SC interface, p-wave pairing can be in- 
duced because inversion symmetry is broken^. Eschrig 
et al. showed that when normal metal is HM, even fre- 
quency pairing would be mostly p-wave^l. Furthermore, 
there are experimental indications of strong proxim ity ef- 
fect between a HM and an s-wave superconductor^^]] 
Here we will show how we can obtain p x + ip y pairing 
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FIG. 1. The heterostructure for obtaining p x +ip y supercon- 
ductivity. The half-metal layer has a 2D electronic structure 
with a single Fermi surface without spin degeneracy. It is 
coupled to the s-wave superconductor through hopping. The 
substrate stabilizes the half-metal crystal structure without 
affecting qualitatively its electronic structure. 

symmetry in a 2D HM when it is coupled to an s-wave 
superconductor only through electron hopping across the 
interface. If the 2D HM has a single Fermi pocket with- 
out spin degeneracy, such p x + ip y pairing will give us 
the TSC with Af = 1. We will show band calculation 
for thin film material that is HM and has a single Fermi 
pocket. We will also discuss the suitable superconductor 
for optimizing this proximity effect and the method we 
can use for detecting the p x + ip y pairing in the HM. 

Basic Model: We consider the model with a s-wave 
superconductor and a 2D half-metal coupled by a weak 
hopping between two systems: 

H = Use + Uhm + H t (1) 

where 

Use =5I( e k ~ Ai') c L c k CT + X]( A k c kt c -l4 + h ' c ')' 

k,CT k 

Uhm =X^( £k ~ M)/kf/kt> 

k 

■Ht =^(tk^/kt c k-+h.c.) (2) 

kcr 

(note k is an in-plane vector). This model provides a 
general mechanism for SC proximity effect. Note that 
while Jit requires k to be conserved in the hopping pro- 
cess, it allows for both spin-flip hopping and dependence 
of hopping on k. 
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FIG. 2. Schematic representation of the HM and SC band for 
weak hopping. Two HM bands, with negative energy portions 
dotted, are due to the artificial doubling of degrees of free- 
dom in the Bogoliubov-de Gennes formalism. In this formal- 
ism, the pairing amplitude and gap require hybridization of 
the 'particle' and 'hole' bands, which can only occur through 
hopping to the SC. 

In this model, the symmetry of Cooper pairs formed on 
the half-metal side, (/-kt/kt)j is determined entirely by 
tk.-fa- I n this model, the only channel for this Cooper pair 
formation is to have the HM electrons with momenta k 
and — k hop to the SC side to form a pair there. This pro- 
cess requires that hopping flips the spin of either one of 
(but not both) the k and — k electrons. Since the Cooper 
pair on the SC side is in the spin-singlet s-wave state, the 
two hopping amplitudes interferes destructively, giving us 



(/-kt/kt) °t ik,tt f -k,n - *-k,tt*k,t4. 



(3) 



with an s-wave multiplicative factor; note the odd spatial 
parity. In the limit of weak hopping, |fk,t<rl ^ l^kd we 
find the pairing amplitude at the HM Fermi surface (i.e. 
for k that satisfies ek — /i = 0) to be22 



(/-kt/kt) 



1 



VW 2 + ICk| 2 



(4) 



where 



Vk — (^k.tt^-k.ti ~~ £k,t^-k,tt)( c -ki c kt)|e 



Ck 



tt 



l*k, 



n\ 



(5) 



as- 



with = - m') 2 + |A k | 2 (|i k ,vl 2 = l*-k, tCT . 

sumed). Eqs. Q and ^ tell us that if the Fermi surfaces 
of the HM and the SC match exactly (i.e. e k — fjf = 
when e k - /i = 0), we have (/-kt/kf) = e l0k (c_ k tc k t) at 
the HM Fermi surface, with e l ^ k being the phase factor of 
ik,tt*-k,t4.-tk,t4/-k,th • Physically, ry k is proportional to 
the amplitude that the k and k HM electrons hop to the 
s-wave SC with opposite spins and form a Cooper pair, 
while ^ is proportional to the amplitude that these elec- 
trons hop to the SC with same spins aligned and therefore 
do not form a Cooper pair. 

Our model gives us not only the pairing amplitude but 
also the pairing gap on the HM side. In the weak hopping 
limit we have been discussing, we obtain 



A 



HM _ ^k _ ^k.tt^-k.t-i- ~ frc.t-l^-k.tt 



At 



(6) 



at the HM Fermi surface^. Note here that |A k | is max- 
imized when we have a perfect Fermi surface matching 
between the HM and the SC. We point out that sinc e our 
HM is 2D, proximity to a SC leads to a pairing gap 2 " 3 24 . 
Eqs.Q and Q tells us this gap A HM comes from the HM 
electron pairing. We therefore conclude that the topolog- 
ical property of SC induced in the HM to be determined 
entirely by A HM , or equivalently, by ik.fo-- 

In the limit of strong hopping, we obtain a much larger 
pairing gap while the pairing symmetry is still deter- 
mined by the hopping as in Eq.([3|. Our strong hopping 
limit requires at the HM Fermi surface — /i = 0, the 
energetics is dominated by the spin-conserving hopping 

i k ,tt) |*k,ffl > E' k and |t k ,ttl > I*k,t4-I- Within this 
model, at the HM Fermi surface (ek — A* = 0), while the 
pairing amplitude remains that of Eq.Q, the pairing gap 
at the HM is now given aP^l 



\HM _ 
^k — 



^k-^k _ *k,tt*-k,t-l- ~~ *k,T4-*-k,TT 



l*k, 



tt 



tt 



(7) 



The HM pairing gap in the strong hopping limit, given 
in Eq.([7]), is much larger than that in the weak hopping 
limit, given in Eq.([6|. This is because in Eq.([7]), A HM 
is proportional to ^J-AtT' wrm e in Eq.([6]), it is propor- 
tional to t^t^/E' 2 . Similar analysis can be applied to 
the proximity of HM with a 3D SC, in which case A^ M 
is still linearly proportional to t^. 

Interface hopping: To show how this HM / SC prox- 
imity effect gives us the p x + ip y pairing, we first discuss 
the symmetry tk,fa-- We note that, since inversion sym- 
metry is broken at the HM/SC interface, th ere w ill be an 
interface Rashba spin-orbit coupling (SOC^P^l 



'H-soc — ha(cr x k) • n<5(n • r). 



(8) 



with n the normal direction of the interface. To account 
for this interface SOC in our model Eq.([2]), we need to 
have a hopping term that has the same symmetry as 
Eq.([8]). Such a hopping term, in the momentum space, 
can be written as 

Ht-soc = tsocy^Fl^ smk y a-cr y sin fc x a)C k + h.c, 



k 



(9) 

where F k = (/kt,/kt) T and C k = (c k t,Ckt) T - However, 
the terms involving /kt can be ignored for the HM/SC 
interface, as these terms involve process occurring at en- 
ergy larger than the minority-spin gap of the HM. 

We can now show explicitly how we obtain the p x + ip y 
pairing from ik,t<r- Assuming the spin-conserving hop- 
ping to be momentum-independent, we obtain 



tk. 



tt 



--to 



b SOC 



[i sin k x a + sin k y a) 



(10) 



We note that for a square lattice in the real space 



representation, the spin-conserving hopping of Eq.(10) 
is 'vertical', i.e. t ^\ fJ*Cif + h.c, while the spin- 
flip hopping is of the 'nearest neighbor' type, i.e. 
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FIG. 3. Crystal structures of two candidate materials, VTe 
and Cr02.The left shows the zinc blende crystal structure of 
bulk VTe. Here, both V and Te forms body-centered cubic 
with the lattice constant of a — 0.6271nm. (The zinc blende 
CrTe has the same structure as the zinc blend VTe, but with 
a slightly smaller lattice constant .) The left shows the ru- 
tile crystal structure of bulk Cr02- The lattice constants are 
a — 0.4421nm and c = 0.2916nm. The distance between the 
nearest O and Cr on the same layer is 0.1817nm. 

t soJ2(ij) ^P[*%]/if c j4. + n - c -> where % gives the chi- 
ral p-wave symmetry to this spin-flip hopping. Inserting 
these tk,t<r into Eq. ^ gives us the chiral p-wave pairing 
on the HM side: 



(/-kf/kf) °t tk,tt*-k,tl — *-k,tt*k.t-l- 

= — 2it t5oc( sm k x a — i sin k y a). 



(11) 



Eq.(10l gives us the Af = 1 TSC in the strong hopping 
limit as well as the weak hopping limilPS] We also note 
that the chiral p-wave is mentioned as the likely pairing 
symmetry if there is intrinsic SC in a HIVpS. 

From the origin of the interface SOC, we can estimate 
of the HM pairing gap to be \A HM \ ~ \A'\(a SO c/W), 
where a$oc is the Rashba SOC of the s-wave SC and 
W is the bandwidth. Physically, when SOC is strong 
for the s-wave SC but weak for the HM, we can expect 
to have the interface SOC. This situation is experimen- 
tally relevant because s-wave SC can exist in materi- 
als with strong SOC, while the spin-polarized ARPES 
indicates complete spin polarization for the HM, as in 
CrOi2£J. Assuming that we have zero effective SOC 
on the HM, we will effectively have the spin-orbit cou- 
pling hopping tsoc ~ (otsoc /W)to induced through sec- 
ond order perturbation. This is sufficient for estimating 
|A ffM |, because, in the strong hopping limit, we find 
A HM ~ (tsoc /h) A> b y inserting Eq.(fil| into Eq.Q. 

We also point out that this HM/SC proximity effect 
can provide us with a means to obtain a multi-domain 
chiral p-wave SC. Eqs.([9]) and (10) show us that if we 
reverse the HM spin polarization, than we will also re- 
verse the chirality of the induced SC. Therefore in a HM, 
a domain boundary between opposite spin polarization 
will also be the domain boundary between the p x + ip y 
and p x — ip y domain when SC is induced. 

Candidate material: We note that the spin-polarized 
(or single-spin) p x + ip y SC with a single Fermi pocket 
gives us the Af = 1 TSCP, but the same can not be said for 
the p x + ip y SC with multiple Fermi pockets^ 1221 . There- 
fore to obtain the Af — 1 TSC, it is important that we 
find a 2D HM with a single Fermi pocket, as we have in 




FIG. 4. (a) The band structure of the VTe-ZnTe(lll) model. 
There is no spin degeneracy in these bands; the red and the 
green dotted curve shows the spin-up and -down bands when 
SOC is absent. The region shaded in orange shows the energy 
range for which we obtain a single Fermi pocket, (b) The sin- 
gle Fermi surface at the energy level O.OeV for ZnTe-VTe(lll) 
model. The green box shows the first Brillouin zone(BZ). (c) 
The band structure of the Cr02 (001) model, done in the 
same way as the VTe-ZnTe(lll) model, (d) The single Fermi 
surface around the F point at the energy level O.OleV for the 
Cr02 (001) model (the green box showing the 1st BZ). 

Eq.{2). That way, the HM / SC proximity effect we dis- 
cussed will give us an equivalent of the single-spin p x +ip y 
SC with a single Fermi pocket 

One candidate material for a 2D HM with a single 
Fermi surface is the zinc blende VTe or CrTe that is two 
atomic layers thick in the (111) direction. We identified 
candidate material through ab-initio band calcula tions 
performed in the frame of density-functional theory^Ml 
with the plane-wave pseudopotential method^. Both 
VTe and CrTe in the zinc blende structure have been 
shown to be half-metallic in band calculation^!! a thin 
film of the zinc blende CrTe has been fabricated in thin 
films by molecular-beam epitaxy^. As we see in Fig. |4j 
two atomic layers of zinc blende VTe (111) on the zinc 
blende ZnTe substrate is a half-metal with a single Fermi 
pocket at the Fermi level; this is due to a mechanisrrp2l 
analogous to th e "po lar catastrophe" in the LaA103 / 
SrTi0 3 interfaccPlMI. As we have a 0.3eV range for the 
Fermi level that gives us a single Fermi pocket, unlike in 
many of the previous proposals for obtaining the Af = 1 
TSO^ we do not require fine-tuning of the HM Fermi 
level. For CrTe, we find a narrower Fermi level range for 
a single Fermi surface (« 0.06eV). 

We also point out that the Cr02 film that is two atomic 
layers thick in the (001) direction comes close to fulfill- 
ing our requirement. Bulk Cr02 has been experimen- 
tally confirmed to be half-metallic 1 26 * 35 ^, and it was also 
shown in experiment to have a strong proximity effect 
to an s-wave SC - NbTiN - that has a relatively high 
T c (~ 14K) and a spin-orbit coupling larger than the 
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SC gap. Since the width of bands near Fermi level for 
NbN is about 7.2e\*™and the atomic spin-orbit coupling 
of Nb is O.leV^, the mechanism discussed here may be 
able to induce a pairing gap up to ~1K. The bulk CrC>2 
crystal, as shown in Fig. [3j has a rutile structure, where 
the Cr atoms form body-centered tetragonal lattice and 
are surrounded by distorted O octahedrsP^. The lattice 
is compressed in the c-axis direction (c/a = 0.66), with 
the lattice constant a — 0.44nm. The band calculation 
for two atomic layer of &O2 (001) in a perfect rutile 
structure gives us a half-metal with one particle pocket 
around the V point and four hole pockets, each with area 
1/4 that of the particle pockets; as shown in Fig. |4j we 
can make hole pockets vanish by raising the energy level 
by less than O.OleV. Note that the bandwidth if very nar- 
row for hole pockets and we obtain a HM with a single 
Fermi surface for a Fermi level range of ~ 0.5eV. 

Detection: Detecting the chiral Majorana edge state 
along with the fully gapped quasiparticle spectrum in the 
bulk of the HM will confirm that we have Af = 1 TSC in 
the HM. The most direct method for detection would be 
measuring with STM the tunneling dl / dV vs. V spectra 
of the HM; indeed this has been used recently to detect 
the SC proximity effect in a If p x + ipy pairing 

is induced, dl/dV is suppressed below V = |A ffM |/e in 



the bulk while along the edge we will see the zero bias 
conductance peak. Another place to look for the zero bias 
conductance peak is the HM domain boundary. This is 
because this will also be the boundary between p x + ip y 
and p x — ip y domain in the SC state and there will be 
two chiral Majorana states running along it. 

In summary, we have shown that we can obtain the 
Af = 1 TSC in a HM through proximity effect with an s- 
wave SC. In a model where the HM is coupled to the SC 
through hopping, the symmetry of SC pairing induced in 
the HM is determined entirely by the hopping term. Due 
to the interface Rashba spin-orbit coupling, the hopping 
term will induce the pairing with p + ip symmetry. In 
order for this p x + ip y pairing to lead to the Af = 1 TSC, 
we need to have a HM with single Fermi surface, and our 
band calculation shows that this can be obtained for a 
very thin Cr02 film. STM measurement can be used to 
verify that TSC is induced in the HM. 
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Appendix A: HM pairing correlation 

Our model Hamiltonian can be solved exactly, as it is 
quadratic: 

H =const. + ^(e k - M)/k t /kf + E ^L^ 

k kcr 

+ XMt{*k,tt( u ik7kt + v ^-ki) 

k 



=*E*lA BdG *k, 
z k 

where * = (/ kt , /! kf , 7kf, 7kL, 7_kf 7- k i) T and 



jjBdG 



Tt 



with /i ffM = r z (e k - 
diag(£ k ,£k, -#k, --Bk), and 



(Al) 



(A2) 



A 



BdG 



"k*k,tt Wk*k,fl 



W k*-k,ti 



W k^-k,tt 



-Vkik.fi 

— Wk^k^tt 



This Hermitian matrix can be 



w k£k,tt 

- u kr k ,u 

(A3) 
diagonalized: 



[/^^ G [/ k = A k , where A k = diag(Ai, . . . A 6 ) and 
[/k is a unitary matrix that becomes an identity matrix 
in the £k,to- — ^ limit. We can also defined the basis 



(/kt, /Ikf 7kt, 7ki, 7l k f 7tki) T that 



*k = c/ k *k 

diagonalize our BdG Hamiltonian. 

We can obtain pairing correlator (/_kt/kt) from this 
diagonalization. Since this diagonalization allows us to 
express the original operators in terms of operators of "J, 

/kt =(E/k)n/kt + (^k)i2./i k t + (^OisTkt 

+ (^k)i47H + (C/k)i57l kf + (Ck)i67li4 

/l kt =(^k) 2 i/kt + (f^Wikt + (^k) 23 7kt 

+ (^k) 2 47ki + (^k)257lkt + (f / k)267lkV 



(A4) 



we can express the pairing amplitude in terms of elements 
of U k : 

(/-kt/kt) =(C/k)l 2 (C/k) 22 + (f/k)l 5 (C/k) 25 + (^k)l 6 (C/k)26 
= (C/k)2l(^k)22 + (tfk)6l(tfk)52 + W(f/k) 62 . 

(A5) 

Since the column vectors of C/k are eigenvectors of H^ dG , 
we can obtain the pairing amplitude by obtaining the 
eigenvectors of H^ dG . In practice, both the weak and 
strong hopping limit enables us to use perturbation the- 
ory for calculating the pairing amplitude. 



1. Weak hopping at the HM Fermi surface 

Here, we assume |ik,t<xl ^ l^kl- At the HM Fermi 
surface, we need to apply the degenerate second order 



perturbation theory since h HM = 0. Therefore, we write 
down eigenvectors only up to terms linear in t ka : 



A« 

-[A( l Wktt + B(4 W- k ,n]/£k 

-[A(*W kn " 5< 4 W-k,tt]/£k 
+[A«) v *t* in - B«<t_ k>n ]/Sk 



(A6) 



(i = 1,2) where AW,fiW are 0(1). In this regime, we 
have 



(/-kt/kt) « [i? (1) ]*S (2) - [A«5«r 



(A7) 



Note that up to quadratic order in t kcr , (A^ l \ B^) T are 
effectively eigenvectors of 

H ef} = ^(-r z Ck - r^Rc^k + rnm % ), (A8) 

^k 



when |t k:tcr | 2 
we need 



Kk,tt 



i-k,t<r| 2 - To justify this approximation 

Mk«k| _ I*k,t4-I 



l*k,ttl ll"kl 



Fk 



> 



|A k | 



(A9) 



at the HM Fermi surface. Within this approximation we 
obtain 



1 



(/-kt/kt) ~ o 



»7k 



2 VlCkl 2 + w 



(A10) 



We also note that the half-metal spectrum is now gapped, 

with the gap of Vl'^+K' iL . However, it is clear that only 
the ?7k part is related to pairing, so we obtain 



(All) 



2. Strong hopping at the HM Fermi surface 

In order to apply perturbation theory in this limit, 
we note that A k /fk,tt ano - *k,t4.Ak,tt are sman variables 
now. It is therefore convenient to rewrite H^ dG in a new 
basis consisting of a spin down electron and a spin down 
hole in the s-wave SC and spin up electrons and holes 
in bonding and anti-bonding states over the HM and the 
s-wave SC. At the HM Fermi surface, this gives us 



fjBdG = U t H BdG u = 



~ h tun (fpaiy 



repair 



A k 



where A k = T z (|u k | 2 - \v k \ 2 )E k , 

h k un = diag(i k ,tt> *-k,tt, _ *k,tt, — *-k,tt) 



, kk| 2 ~K| 2 F 

H ^ -C/k 



1 -1 

-1 1 
-1 1 

1 -1 



, (A12) 
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fl air = y/2 



and 



*k,n/ 2 <«k£k t* Kn /2 -ulv k E u 
UkV^E k -t_k,tj./2 -UkV^Ek -t-k M /2 

(A13) 



s/2 



v'2 



1 

x/2 

i'k 
V2 
_^k 
n/2 



1 



"k 



V2 



1 

n/2 



\/2 
v/2 



«k 



(A14) 



(note that chose a gauge so that ik,tt ^ s rea ^ positive). 

Our strong-hopping mode in 2D1 gives us three spin- 
polarized bands. We can see this from examining the 
column vectors of U. Note that the first and third 
columns represents the electron in the spin-up bond- 
ing and anti-bonding bands respectively, the second and 
fourth columns represents the hole in the spin-up bond- 
ing and anti-bonding bands respectively, and the last two 
columns represents the electron and the hole in the spin- 
down band. For ik,tt = ^-k,tt' we nave three pairs of a 
electron state and a hole state degenerate in energy - two 
pairs for spin-up and one pair for spin-down. Therefore, 
we can again apply degenerate second order perturbation 
theory with three sets of bases 



AW 



Kl 2 -I"kl 2 



AW 



4 *k,TT 

Kl 2 -M 2 -Ek 

4 £-k,TT 



V'2 * 



B k 



t— k,tt 
1 t-kjti 



y/2 t- 



(A15) 



4 *k,tt 

_ D (i) |n k | 2 -|fk| 2 E k 

4 t_k,fr 



B k 
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(A16) 



K 



W 



\/2 *k,ft 



B k 



- y(« 



tt 

k,T4 



v'2 t-k.ff 

B k 



-X(V2u k <dk, tt 

'^^t-^^k^ 

) J_ *-k,T4. 
v'2 *-k,TT 



AW 
y« 



(A17) 



We can show from the effective Hamiltonians for the 
above three bases that the chirality of the spin-down pairs 
will be opposite to that of the two sets of spin-up pairs. 



The effective Hamiltonians of and can be writ- 
ten 



£k 



(r 2 Ck + r^Re^k - r^In^k) 



i l*k,ft 
- u\ 

while for , the effective Hamiltonian is 



1 |A k | 2 i|^,nl 2 , 1(4-/0 

2 |i k ,nl 2 2 |i k , n l 2 



/\21 



Knl 2 

(A18) 



2S k 



^(Ak-T^k-T^k). (A19) 



Because the tunneling process required for pairing the 
spin-down pairs is the reverse of what is required for pair- 
ing the spin-up pairs, the chirality of the spin-down pairs 
is opposite of that of the spin- up pairs. Eqs.(A18) and 



( A19 1 show that the HM / SC heterostructure overall has 
the N = 1 TSC. We believe that this result should hold 
even in the case the SC is 3D. 

These effective Hamiltonian tells us that we have giv- 
ing us 



A 



HM 



>/k 



|<k, 



(A20) 



since only the spin-up pairs involve the HM. 



Appendix B: Polar catastrophe in zinc-blende 
VTe/ZnTe (111) 

In this Appendix, we discuss the physical origin of the 
"polar catastrophe" in the zinc-blende VTe/ZnTe (111) 
that gives us a single Fermi surface at O.OeV. To explain 
this, we first need to clarify what crystal structure we 
used for the band structure calculation. To maintain 
the zinc blende structure of VTe, we always need a zinc 
blende substrate with a good matching lattice constant. 
The zinc blende ZnTe not only has the lattice constant 
matching within 1% but also has been used experimen- 
tally in growing the zinc blende CrTe. For VTe, growing 
either of them on the (111) side of the zinc blende ZnTe 
is sufficient in obtaining both half-metallicity with a sin- 
gle Fermi surface and the two atomic layer zinc blende 
(111) structure with very little structural distortion. The 
crystal structure of this heterostructure can be a series 
of monatomic triangular lattice in the sequence V-Te- 
Zn-Te-Zn-Te-... with all Zn atoms sitting at the center 
of tetrahedron formed by Te atoms while Te atoms adja- 
cent to V is one short of forming tetrahedron. (The zinc- 
blende CrTe/ZnTe (111) has essentially the same crystal 
structure.) 

We now explain how we effectively have three and one- 
half electrons per unit cell of the VTe at the top. This 
is crucial because having a fractional number of electron 
per unit cell is necessary for a half-metal to have a single 
Fermi pocket at the Fermi level. We first note that ZnTe 



G 



is an insulator (with a gap of 2.3 eV) due to covalent 
bonding between Zn and Te takes up all valence elec- 
trons of this atom, with Zn contributing 2 electrons per 
atom, or equivalently, 2/4 = 1/2 electrons per bonding 
with Te. We need to note here that the electron configu- 
ration of a V atom differ from that of a Zn atom due to 
addition of three 3c? electrons. Since each V atom is miss- 
ing one Te atom, we get 1/2 electrons for every dangling 
bond from V. While these electrons from dangling bonds 
cannot flow into the bulk ZnTe, the VTe film is metallic 
due to V 3d electrons. In short, the dangling bonds and 
3c? electrons conspire to give us a half-metal with a sin- 
gle Fermi pockets, and it is from these four orbitals that 
we get one partially filled bands and three filled bands 
close to the Fermi level in Fig 4 (a). As we have pointed 
out, this is analogous to the "polar catastrophe" at the 
LaA10 3 / SrTi0 3 interface. 

We have set up our calculation so that this "polar 
catastrophe" would not lead to any breakdown. In our 
calculation, when we used a ZnTe of finite thickness, the 
dipole moment between Zn and Te atoms give rises to 
an electric field in our heterostructure. Since this effect 
decreases when we increased the thickness of the ZnTe 



layer, we expect the effect to go away for a sufficiently 
thick ZnTe layer. For our band calculation, we added 
an extra atomic layer of the H atoms - which gives a 
VTe/ZnTe/H heterostructure - solely to cancel out this 
artificial electric field. In Fig 4 (a) we have omitted the H 
band, which do not mix with any VTe band due to both 
the insulating gap and the thickness of the ZnTe layer. 



Appendix C: The method for ab-initio calculations 

In this work, all the ab-initio calculations are per- 
formed by the BSTATE(Beijing Simulation Tool of 
Atomic Technology) package with the plane-wave pseu- 
dopotential method. The Perdew-Burke-Ernzerhof type^ 
generalized gradient approximation(GGA) is used for 
exchange-correlation energy. Due to 2D models, the k- 
point grid is taken as 12 x 12 x 1 for general self-consist cal- 
culations, and the kinetic energy cutoff is fixed to 340eV. 
All the atoms on the surface part are fully relaxed for our 
free-standing slab zb-VTe/ZnTe (111) model. However, 
for our Cr02 model, we assumed that the perfect rutile 
structure and applied no relaxation. 
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